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ABSTRACT 

We study the four double neutron star systems found in the Galactic disk in 
terms of the orbital characteristics of their immediate progenitors and the natal 
kicks imparted to neutron stars. Analysis of the effect of the second supernova 
explosion on the orbital dynamics, combined with recent results from simulations 
of rapid accretion onto neutron stars lead us to conclude that the observed 
systems could not have been formed had the explosion been symmetric. Their 
formation becomes possible if kicks are imparted to the radio-pulsar companions 
at birth. We identify the constraints imposed on the immediate progenitors of 
the observed double neutron stars and calculate the ranges within which their 
binary characteristics (orbital separations and masses of the exploding stars) 
are restricted. We also study the dependence of these limits on the magnitude 
of the kick velocity and the time elapsed since the second explosion. For each 
of the double neutron stars, we derive a minimum kick magnitude required for 
their formation, and for the two systems in close orbits 10 R©), we hnd it 
to exceed 200km/s. Lower limits are also set to the center-of-mass velocities of 
double neutron stars, and we hnd them to be consistent with the current proper 
motion observations. 


Subject headings: stars: neutron - pulsars: general - supernovae: general 
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1. Introduction 

The Double Neutron Star (DNS) system containing PSR 1913+16, discovered in 
1975 by Hulse & Taylor, was the hrst binary pulsar for which observational data strongly 
suggested that both components were compact objects. Following PSR 1913+16[|, additional 
DNS candidate systems have been discovered: PSR 2303+46 (Stokes, Taylor, & Dewey 
1985), PSR 2127+llC (Anderson et ah 1990), PSR 1534+12 (Wolszczan 1991), and PSR 
1518+4904 (Nice, Sayer, & Taylor 1996). The high orbital eccentricities of these systems 
separate them from all other binary radio pulsars, which are found in circular orbits (e.g., 
van den Heuvel 1995). Measured and derived orbital and pulsar properties of these hve 
systems are given in Table 1. 

One of the two members in each binary is a neutron star since it is detected as a 
radio pulsar. Using a variety of tests, the nature of the companion to the pulsar has been 
investigated in detail for PSR 1913+16 (Taylor & Weisberg 1982, 1989) and to a lesser 
extent for the other systems. Extended companions are excluded based on their radial size, 
their implied periastron advances, and their apparent visual magnitudes, while helium stars 
of the same mass are even brighter (Crane, Nelson, & Tyson 1979; van Kerkwijk 1996; 
Webbink 1996). For three of the systems, white dwarf companions of various temperatures 
can potentially be ruled out using deep optical observations (see Table 2). For the two most 
distant systems, PSR 1913+16 and 2127+llC, there is no current observational test that 
can rule out white dwarf companions. However, the evolutionary history and high orbital 
eccentricities of these systems require that the pulsar companions were formed in supernova 
explosions that occurred after the pulsar spin-up phase (Srinivasan & van den Heuvel 1982). 
This provides indirect evidence that all of the pulsar companions are either neutron stars 

^Hereafter, we will refer to the double neutron star systems using the designation of their 


radio pulsar component. 
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or black holes. 

The evolution of double neutron stars after the second supernova explosion is governed 
only by angular momentum losses due to emission of gravitational waves, and is not 
complicated by effects of binary-star evolution (e.g., nuclear evolution, mass transfer, stellar 
winds). These systems, therefore, constitute an ideal diagnostic for neutron star kicks. 
Natal kicks are thought to be imparted to neutron stars at the time of their formation as 
a result of some asymmetry in the supernova explosion. Evidence for the existence of such 
kicks originates in observations of radio pulsar velocities (Caraveo 1993; Harrison, Lyne, 

& Anderson 1993; Frail, Goss, & Whiteoak 1994; Lyne & Lorimer 1994), 0,B runaway 
stars (Stone 1991), and X-ray binaries (van den Heuvel & Rappaport 1987; Kaspi et ah 
1996), and estimates of their average magnitudes vary from ~ 100 to ~ 500km/s (Dewey & 
Cordes 1987; Bailes 1989; Fryer, Burrows, & Benz 1997). 

In this paper, we thoroughly investigate the four double neutron star systems found 
in the Galactic held to identify the allowed binary characteristics of their immediate 
progenitors and to constrain the nature of the natal kick required for their formation. In § 2, 
we examine the possibility that DNSs are formed via symmetric supernova explosions. Based 
on arguments related to orbital dynamics and, in contrast to previous work, accounting for 
orbital evolution due to gravitational radiation, we hnd that the hrst neutron star in the 
observed systems always lies within the envelopes of their exploding companions. In the 
past, similar conclusions were reached for some of the systems (Burrows & Woosley 1986; 
Yamaoka, Shigeyama, & Nomoto 1993), although the fate of the engulfed neutron stars 
was not addressed in any detail. We use recent results on neutron-star common-envelope 
phases, and conclude that the second supernova explosions in DNS systems could not have 
been symmetric. In § 3, we account for the effects of kicks imparted to nascent neutron stars 
and hnd that the observed systems can then be reproduced. Previous studies (Burrows 
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& Woosley 1986; Wijers, van Paradijs, & van den Heuvel 1992; Yamaoka, Shigeyama, & 
Nomoto 1993) have considered kicks of specific magnitude or direction or have indicated 
plausible values of progenitor parameters and kick magnitudes. Here, we are able to (i) 
identify the ranges of values within which the properties of the immediate progenitors of the 
observed DNSs are restricted, and (ii) derive an absolute lower limit to the magnitude of the 
kick velocity needed for their formation. We conclude with a discussion of the implications 
of our results for the intrinsic kick velocity distribution and the center-of-mass velocities of 
double neutron stars. 


2. Symmetric Explosions 

The existence of double neutron stars has instigated studies of their origin and a 
variety of evolutionary sequences have been proposed for their formation (Flannery & 
van den Heuvel 1975; Webbink 1975; Smarr & Blandford 1976; Srinivasan & van den 
Heuvel 1982; van den Heuvel, Kaper, & Ruymaekers 1994; Brown 1995; Terman & Taam 
1995). Regardless of the details of prior evolution, all mechanisms converge to the same 
conhguration before the formation of the second neutron star: a binary consisting of a 
neutron star and a helium star in a circular orbit. The helium star then collapses to form 
the companion to the observed radio pulsar. A possible exception is PSR 2127+llC, which 
is found in the dense stellar environment of a globular cluster and may have formed via 
dynamical interactions (capture or exchange: Anderson et al. 1990). Therefore, we restrict 
our analysis to the four DNS systems found in the Galactic held. 

In this section, we assume that the supernova explosion of the helium star is symmetric, 
its only effect being an instantaneous mass loss from the system. Since the companion of the 
exploding star is a neutron star, we may neglect any interaction between the radio pulsar 
and the supernova ejecta. Combining the results for the pre-SN conhguration with recent 



developments on common envelope evolution of inspiraling neutron stars, we conclude that 
the four double neutron star systems found in the Galactic held could not have been formed 
if supernova explosions were symmetric. 


2.1. Orbital Dynamics 

In order to obtain the post-SN parameters for each of the known systems, we must 
account for orbital evolution due to gravitational wave emission (Wijers, van Paradijs, & 
van den Heuvel 1992 recognized the distinction between present and post-SN characteristics, 
but did not include this effect in their analysis). For each system, the extent of orbital 
evolution depends on the time, Tsn, elapsed since the supernova explosion. Since the radio 
pulsar has been recycled prior to the supernova explosion, its age sets an upper to the 
value of Tsn- The difference between the two times is of the order of the lifetime of the 
helium star (~ 10® yr; see Habets 1985). The “characteristic age” of the pulsar dehned 
as Tc = Ps/2Ps (where Pg and Pg are the pulsar spin period and its derivative) is often 
used as a measure of the pulsar age. However, this estimate is strongly model dependent 
(assumptions about the form of the braking law, the pulsar magnetic held and the initial 
spin period), and hence the upper limit on Tsn, may be considerably different than the 
value of Tc (e.g., Lyne & Graham-Smith 1990). Therefore, in our calculations, we treat Tsn 
as a free parameter ranging from O.Itc to IOtc. We use the expressions derived by Junker & 
Schaefer (1992) for orbital evolution due to gravitational radiation to integrate backwards 
in time, and calculate the post-SN orbital separations. A, and eccentricities, e, from the 
observed values (Hcur,ecur)- This effect is strongest for the two systems with short orbital 
periods, PSR 1913-1-16 and PSR 1534-1-12 (see Figure 1). 

For a given set of post-SN characteristics, we can calculate the corresponding 
parameters just prior to the explosion. We denote these parameters with a subscript “o”. 



The radio pulsar mass, Mp, remains the same before and after the explosion. 


2.1.1. Circular Pre-SN Orbits 

In the case of circular pre-SN orbits, the orbital separation, Ao, and mass. Mo, of the 
exploding star are given by (see also Blaauw 1961; Boersma 1961): 


Ao 

Mp + 2 Mo - Mo 

(1) 

A 

~ Mp + Mo 

Mo 

= 6 (iWp -/Wc) Mq. 

(2) 


According to the current formation mechanisms, the exploding star is an helium star. 
During the course of their evolution, helium stars less massive than ~3.5M0 expand 
signihcantly. We use a fit for the maximum helium-star radius, i?He,max, as a function of its 
mass, Mne, (Kalogera & Webbink 1996) based on helium star models calculated by Habets 
(1985) (which are consistent with the models of Woosley, Langer, & Weaver 1995): 

\ 3.0965 - 2.013 log Mne Mue < 2.5 Mq 

logi?He,max = < (3) 

[ 0.0557 (logMne - 0.172)-2-5 M^e > 2.5 Mq. 

Using equations (1), (2), and (3) we calculate the pre-SN separations, Ao, progenitor 
masses of the pulsar companions. Mo, and the ratios of orbital separations over helium-star 
maximum radii, Ao/Rue,max, as functions of Tsn. The results for PSR 1913-1-16 and PSR 
1534-1-12 are shown in Figure 2, while the ranges of values for all systems are given in Table 
3. It is evident that for all observed DNSs and for a wide range of Tsn, the stellar radius 
of the helium star exceeds the orbital separation. In other words, if the explosion were 
symmetric, we hnd that in the progenitors of all the observed systems, the radio pulsar 
would he within the envelope of its companion at the time of the supernova explosion (see 
also Burrows & Woosley 1986; Yamaoka et ah 1993). 



2.1.2. Eccentric Pre-SN Orbits 


All current formation scenarios for double neutron star systems lead to circularization 
of the pre-SN orbits as the hrst neutron star gets recycled. However, the pulsars may not be 
recycled or there may be a way to recycle the pulsar without circularizing the orbit of the 
system. Hence, we examine our conclusion that the pulsar would be in a common envelope 
phase at the time of the explosion, in the case of eccentric pre-SN orbits. 


In this general case, the relations between the pre- and post-SN parameters have been 
derived by Hills (1983): 


A _ {Mo + Mp) - {Mo - Mo) 

A„ {Mo + Mp) -2 {Ao/r) {Mo - ’ 


and 


, , / Ao Mo- Mo 

e = {l-{l-el) 1 - 2^ 


r Mo + M^ 


1 - 


Mo-Mo 

Mo + Mp 


-2\ 1/2 


( 4 ) 


( 5 ) 


where r is the distance between the two stars with masses Mo and Mp, at the time of the 
supernova explosion. There are four unknown pre-SN parameters: the orbital separation, 
Ao, the helium-star mass. Mo, the pre-SN eccentricity, Co, and the distance, r, between the 
two stars at the time of the explosion. Since we only have two equations, we treat two of the 
above quantities, Co and Mo, as parameters and solve for Ao and r. Note that r can acquire 
values only within a restricted range, between periastron and apastron of the pre-SN orbit. 


The ratio of the periastron distance to the maximum radius for the helium star as a 
function of Mo and for a complete range of pre-SN eccentricities is shown in Figure 3 for all 
four of the observed systems (see also Table 3). We see that although for PSR 2303-1-46 and 
1518-1-4904 the ratio Ao(l — eo)/i?He,max might exceed unity, for PSR 1913-1-16 and 1534-1-12 
the pulsar would he well within the envelope of the helium star. Therefore, if the supernova 
explosions were symmetric, the pulsars in at least two of the observed DNS systems would 
experience common envelope evolution, even in the extreme case of eccentric pre-SN orbits. 



2.2. Neutron Stars in Common-Envelope Phases 


For a symmetric explosion, the immediate progenitors of double neutron star systems 
must evolve through a common-envelope phase. In this section, we use recent results on 
rapid accretion onto neutron stars to examine their fate in the envelopes of evolved helium 
stars. We conclude that survival of neutron stars through such phases is not possible. 

Accretion of matter onto the neutron star in common envelope phases has been 
recently studied by various investigators (Chevalier 1993, 1996; Brown 1995; Fryer, Benz, & 
Herant 1996). In common envelope situations, the Bondi-Hoyle accretion rate (Bondi 1952) 
typically exceeds ~ lO“^M0yr“^, which is much greater than the photon Eddington limit 
for neutron stars. For such high accretion rates, neutrino emission dominates the cooling 
processes. Accretion is not limited by the photon Eddington limit, but instead by the 
corresponding limit for neutrinos, which is 20 orders of magnitude greater than the photon 
Eddington limit. For the calculated Bondi-Hoyle rates, neutrino emission carries away the 
energy released by accretion, a steady state accretion flow at high rates is achieved, and the 
pulsar collapses into a black hole on time scales of hours or months. 

Depending on the conditions in the common envelope, other physical processes can 
facilitate the survival of the neutron star. Chevalier (1996) showed that the angular 
momentum of the accreted material may prevent rapid accretion and hence collapse of the 
neutron star into a black hole. Also, Fryer et al. (1996) found that the neutrinos emitted 
close to the neutron star may be able to heat a region of the infalling atmosphere, and 
ultimately drive a small scale explosion which temporarily halts mass accretion. A series 
of such explosions may also prohibit the collapse of the neutron star into a black hole if 
a sizable amount of the atmosphere can be blown off of the neutron star instead of being 
accreted. Both of these effects depend sensitively on the Bondi-Hoyle infall rate onto the 
neutron star and they become more important as the infall rate decreases. The envelopes 
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of hydrogen-rich giant stars (with masses ~ 10 Mq) are less dense (by at least one order of 
magnitude) than those of helium-rich giants and have correspondingly lower infall rates. 
Although, in helium giants, the infall rates are too high for angular momentum to limit the 
flow, it may well limit the accretion for some hydrogen giants (see Table 1 in Fryer et ah 
1996). 

In this section, we calculate the time, Tcou, needed for a neutron star to collapse into a 
black hole (for zero and maximum angular momentum of the material), and compare it to 
the evolutionary time, Tgv, from maximum radial extent of the helium star to its explosion. 
Only if Tev is shorter than Tcoii the helium star will explode in a supernova before the 
neutron star collapses into a black hole. 

For each of the observed systems and under the assumption of symmetric supernova 
explosions, we have derived the mass of the progenitor of the second neutron star as well 
as the pre-SN orbital separation, that is the position of the pulsar in the envelope of its 
companion at the time the core of the helium star collapses. Using helium star models 
by Woosley (1996) we calculate the Bondi-Hoyle accretion rates, Mb, appropriate for the 
physical conditions in the helium envelopes (see Table 4). We have assumed that the 
neutron star velocity in the envelope is equal to the pre-SN relative orbital velocity. These 
Bondi-Hoyle rates are too high for the explosions found by Fryer et ah (1996) to occur. In 
the extreme case that angular momentum has no effect on accretion, matter is accreted by 
the neutron star at the Bondi-Hoyle rate. Thus, we can calculate the corresponding time 
scale for collapse, Tgou, into a black hole, assuming that the amount of mass required for 
collapse is ~ 1 Mq (see Table 4). 

Matter accreted onto the neutron star is expected to have angular momentum due to 
density and velocity gradients in the envelope. As a result, infall of material is halted at all 
polar angles except in a cone around the poles dehned by the angular momentum symmetry 
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axis. For the extreme case of negligible angular momentum transport, we derive a lower 
bound to the accretion rate onto the neutron star. For a given value of the specihc angular 
momentum, j, there is a critical polar angle, 6c, within which angular momentum does not 
affect the inflow appreciably. This angle is set by the balance between the centrifugal and 
gravitational forces at the neutron star radius (in the Newtonian limit): 

fsin^dc GMns 

-'^NS -'^NS 

where and i?NS ai'e the gravitational mass and radius of the neutron star, assumed 
here to be equal to 1.4 Mq and 10® cm, respectively. For the specihc angular momentum we 
use the analytical estimates by Ruffert & Anzer (1995): 

j = ^{Q(^v - (^p)VR b, (7) 

where V is the velocity of the neutron star, Rb is the Bondi-Hoyle accretion radius (Bondi 
1952), and e^,, are the velocity and density gradient parameters, as dehned in Ruffert & 
Anzer (1995). We note that their numerical calculations indicate that the above expression 
(eq. [7]) overestimates the amount of angular momentum accreted onto the neutron star, 
and therefore the corresponding mass accretion rate we derive is an underestimate. Using 
the helium-star models by Woosley (1996) we calculate j and 6c- For accretion restricted 
within the polar regions, we obtain an extreme lower limit to the accretion rate onto the 
neutron star and a corresponding upper limit to the characteristic time scale for collapse, 

Fcoll- 

The neutron star will avoid the collapse into a black hole only if Tcou is longer than 
the evolutionary time interval, Tev, between the time the helium star acquires its maximum 
radial extent and the time it explodes as a supernova. For all of the observed systems 
(with the possible exception of PSR 1518-1-4904), the collapse time scale is many orders of 
magnitude shorter than the evolution time scale, which clearly indicates that the pulsar 
becomes a black hole long before the SN explosion. 
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In view of the above analysis, we conclude that the pulsar cannot survive a common- 
envelope phase with its helium star companion. Nevertheless, such a phase could not be 
avoided if the explosion were symmetric, and hence would not lead to the formation of the 
double neutron star systems. 


3. Asymmetric Explosions 

We have already shown that if supernovae explosions were symmetric, the hrst-born 
neutron stars would have collapsed into black holes during a common-envelope phase with 
their helium-star companions, thus aborting DNS formation. In this section, we hnd that 
asymmetric explosions, during which kicks are imparted to nascent neutron stars, allow 
the formation of the observed double neutron star systems (see also Yamaoka et ah 1993). 
The pre-SN parameters of their progenitors are still restricted within relatively narrow 
ranges, which depend on the magnitude of the kick velocity, the time elapsed since the SN 
explosion, and the eccentricity of the pre-SN orbits. For each of the systems we are able to 
evaluate the minimum kick magnitude required for its formation. 


3.1. Circular Pre-SN Orbits 


First, we study the case of circular pre-SN orbits, consistent with all current DNS 
formation mechanisms. Using the orbital energy and angular momentum equations for 
eccentric orbits, and for a kick velocity of specihc magnitude and direction, expressions for 
the pre-SN orbital separation, A^, and mass of the helium star, Mq, as a function of the 
post-SN binary characteristics can be derived (e.g.. Hills 1983; Kalogera 1996): 


A 


+ 2VfcU,cos0 = G(Mp + M,) A - , 

UfcW^cosV + (f4cos0 + K)'l = G(Mp + Me)A(l-e2), 


( 8 ) 

(9) 
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where Vr is the relative orbital velocity before the explosion, 


K = 


G(M, + Mp 

An 


1/2 


( 10 ) 


14 is the magnitude of the kick, 6 is the polar angle of the kick with respect to the 
pre-SN orbital velocity of the helium star relative to the pulsar, and 0 is the corresponding 
azimuthal angle, so that 0 = 0 represents a plane perpendicular to Aq. For each of the 
observed DNS systems, the post-SN parameters are known. First, we calculate the orbital 
evolution due to gravitational radiation assuming Tsn = Tc (Table 1), and discuss how 
variations in Tsn affect our results at the end of this section. For a specihed kick magnitude 
in any direction, we calculate the two unknown pre-SN parameters, Ao and Mq. 


The region of allowed pairs of values is restricted by a set of constraints that DNS 
progenitors must satisfy. The first of the constraints arises from the fact that the post-SN 
orbit must include the position of the two stars at the time of the explosion (e.g., Flannery 
& van den Heuvel 1975; Yamaoka et ah 1993; Kalogera 1996), hence the pre-SN separation 
must be within the bounds set by the periastron and apastron distances in the post-SN 
orbit 

A{1-e) < Ao< A{l + e). (11) 

These two limits follow from equations (8) and (9) with the constraint that cos^ 0 < 1, 
where the equality corresponds to kicks restricted to the plane perpendicular to Aq. 


A second constraint is set by demanding that cos^ 0 remains always positive. Using 
equations (8) and (9) we hud that for a specihc pre-SN separation the mass of the helium 
star is restricted to: 


Mn < 


2e{Mp + M^ 


2a [2 0:^(1 — e^) — k] — 4 q; 2 (i _ e2)V2 [q-s _ g 2 ^ _ /.j 


1/2 


( 12 ) 


where 


V^A 


G{Mp + 


+ 1 


k = 2a — 


(13) 
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and a = AjA^. The extreme case of cos^ 0 = 0 corresponds to kicks restricted to the plane 
of the orbit before the explosion. 

The hnal constraint arises from the condition that the pulsar must avoid a common 
envelope phase with the helium star. The ranges of allowed values for A^ and are then 
restricted by the condition that the helium star must £t within its Roche lobe when its 
radius is at a maximum (eq. [3]). At a given value of Aq, a lower limit is derived for the 
mass of the helium star, since the maximum radii decrease with increasing helium-star 
mass. 

The limits imposed on pre-SN orbital separations, Aq, and helium-star masses. Mo, are 
shown in Figures 4 and 5 for the two observed systems for which the individual masses are 
accurately known (PSR 1913-1-16 and 1534-1-12). For the close DNS systems, the value of Ao 
is restricted in the range ~ 3 — 4.5 R© and that of Mo must exceed ~ 4 M©, whereas, for the 
two wide systems (PSR 2303-1-46 and 1518-1-4904), it is Ao ~ 20 — 50R© and Mo ^ 3Mq. 
The solid thick vertical line corresponds to the geometrical constraint (eq. [11]) and sets an 
upper limit on the pre-SN orbital separation. Although there also exists a lower limit to 
Ao, its value is too low to restrict the parameter space. The thin lines correspond to kicks 
imparted in the pre-SN orbital plane (cos^ 0 = 0), and thus to maximum helium star masses 
(eq. [12]); of the three limiting constraints this is the only one that depends on the assumed 
kick magnitude. As the magnitude of the kick velocity increases, the relative orbital velocity 
of allowed DNS progenitors also increases, and systems with shorter separations and more 
massive helium stars are included in the progenitor parameter space (see Figure 4). For the 
two close systems, PSR 1913-1-16 and 1534-1-12, and for a DNS age equal to Tc, the pulsar 
companions must have been imparted a kick of at least 260 and 220km/s, respectively. For 
the wider systems, PSR 2303-1-46 and 1518-1-4904, the minimum kick magnitudes are lower, 
70 and 50km/s, respectively. 
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The parameter space allowed to DNS progenitors depends also on the time, Tsn, 
elapsed since the supernova explosion, during which signihcant orbital evolution due to 
gravitational wave emission possibly occurs. The limits on the ranges of pre-SN parameters 
are altered with the exception of the limiting curve related to the size of evolved helium stars 
(thick solid lines in Figures 4 and 5). As the value of Tsn increases, post-SN separations 
and eccentricities increase too, and hence the maximum allowed pre-SN separations are 
shifted to higher values (eq. [11]). The consequent decrease in relative orbital velocity 
before the explosion also results in a decrease of the minimum kick magnitude required 
for the formation of the observed systems. For Tsn = 10 Tc (or Tsn = 0) the minimum 
kick magnitude required decreases to 170km/s and 190km/s (increases to 290km/s and 
225km/s) for PSR 1913-1-16 and 1534-1-12, respectively. For the two wide systems (PSR 
2303-1-46 and 1518-1-4904) the effect of orbital evolution due to gravitational wave emission 
is negligible, and therefore the minimum kick magnitudes remain unaffected. 


3.2. Eccentric Orbits 

Although none of the current formation mechanisms requires an eccentric orbit for 
the pulsar-helium star binary, as in the case of symmetric explosions, we explore such a 
possibility and its effect on the pre-SN binary parameters. In the case of eccentric orbits, 
the pre-SN parameter space becomes four-dimensional: the eccentricity, Co, and the position 
along the orbits at the time of the explosion expressed by the eccentric anomaly, Eq, are 
added to the orbital separation, Aq, and helium-star mass, Mq. We sample the complete 
range of values for both Cq (0 to 1) and Eg (0 to tt), and calculate the outer envelope of the 
limits imposed on pre-SN separation and mass of the exploding star, based on the three 
constraints discussed above. These limits are shown in Figure 6 for PSR 1913-1-16 and 
PSR 1534-1-12 for two different values of the kick magnitude and for post-SN parameters 
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corresponding to Tsn = Tc- The area of the parameter space allowed to DNS progenitors is 
enclosed by the limiting curves. The lower limit on represents the constraint that the 
helium star at its maximum extent must £t within the orbit. The upper limit corresponds 
to kicks restricted to the orbital plane (cos^^ = 0). The condition cos^^ = 1 is no longer 
restrictive for the progenitor parameters. Since the explosion is instantaneous, the distance, 
r, between the two binary members before and after the supernova remains unchanged: 

r = Ao{l — CoCosEo) = A{l — ecosE). (14) 

In the limit of nearly parabolic pre-SN orbits and for = 0, no upper limit can be set 
to the orbital separation, A^. As the magnitude of the kick increases, wider and more 
eccentric progenitor orbits and more massive helium stars are allowed. However, we note 
that helium stars are known to lose mass in strong stellar winds and thus their masses at 
collapse are expected to be relatively small SMq) (Woosley, hanger, & Weaver 1995). 
Despite the expansion of the parameter space allowed to DNS progenitors, the minimum 
kick magnitude required to form the observed systems remains practically unaltered. 


4. Discussion 

A careful account of the change in orbital characteristics of binary systems experiencing 
symmetric supernova explosions, of the sizes of helium stars approaching collapse into a 
neutron star, and of recent results concerning the fate of neutron stars within the envelopes 
of helium stars, led us to conclude that the observed double neutron star systems in the 
galactic disk could not have been formed if the SN explosions forming the pulsar companions 
were symmetric. Natal kicks imparted to neutron stars are required to explain the observed 
parameters of DNS systems, even if they have experienced signihcant orbital evolution 
after their formation due to gravitational wave emission. In fact, only when kicks have 
magnitudes exceeding some minimum value can the progenitor orbits be sufficiently wide 
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to accommodate evolved helium stars, and yet produce the small separations measured in 
these systems. For the observed double neutron stars in close orbits, the minimum kick 
velocities required are in excess of 200km/s. 

The analysis presented here has enabled us to constrain the orbital parameters of the 
immediate progenitors of the observed DNS systems. However, obtaining the most probable 
of the parameters within the allowed ranges requires modeling of the evolutionary history 
of the DNS progenitors up to the formation of the pulsar/helium-star binaries. We have 
used the population synthesis code (Monte Carlo) described in Fryer et al. (1997) to obtain 
an estimate of the most probable kick velocities, and hence progenitor parameters, for the 
observed systems. The standard DNS formation mechanism (Srinivasan & van den Heuvel 
1982) is modeled, including the constraint that the helium stars must be accommodated in 
their orbits prior to collapse. The formation rate of DNS systems with separations ^ SRq 
peaks at 200km/s (5 x 10“®yr“^), and falls off sharply beyond 300km/s. Although higher 
kicks are allowed by our analysis, they disrupt most of the binaries and do not contribute 
signihcantly to DNS formation. For double neutron star systems with separations ^ 50 R©, 
the formation rate (~ 10“®yr“^) peaks at kicks lower than lOOkm/s, and is consistent with 
the rate inferred observationally (van den Heuvel 1995). These results indicate that kick 
velocities close to the derived lower limits are favored. Conseqnently, progenitor orbital 
separations and masses of ~ 4.5 Rq and ~ 4.5 Mq are favored for the close systems (PSR 
1913+16 and 1534+12). Separations of ~ 50R© and ~ 30R© are more probable for the 
progenitors of PSR 2303+46 and PSR 1518+4904, respectively, with helium-star masses of 
r'U 3 Mq in each case. 

We use our knowledge of the parameters of DNS progenitors to calculate the 
center-of-mass velocities, Vcm, acquired by the systems after the explosion (equation 34 in 
Kalogera 1996), for a specihc kick magnitude. For each observed system, the lowest possible 
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value of VcM corresponds to the minimum kick magnitude needed for its formation, and 
to the progenitor binary with the lowest orbital relative velocity (widest orbit and least 
massive helium star). For the two systems in close orbits (PSR 1913+16 and 1534+12) 
we hnd minimum center-of-mass velocities of 200km/s and 225km/s, respectively. Since 
the results of our synthesis calculations of DNS formation favor low kick magnitudes, 
these minimum center-of-mass are also the most probable. We note that if the age of the 
system is longer than Tc, for example is equal to 10 Tc, then the center-of-mass velocities are 
reduced to lOOkm/s and 185km/s, respectively. For the two wider systems, the minimum 
values of Vcm are considerably lower, Vcm ~ 50km/s. Measurements of proper motions of 
the observed DNSs have been reported in the literature the past few years (Taylor 1993; 
Arzoumanian 1995; Frail 1996; Nice 1996). The inferred values of transverse velocities for 
the two close systems, PSR 1913+16 and 1534+12, are f«110km/s and ?«160km/s. For 
PSR 1518+4904 only an upper limit of 40km/s can be set, while for PSR 2303+46 there is 
no signihcant measurement. We note, however, that these values depend sensitively on the 
distance estimates. Moreover, these measurements include only the transverse component 
of the space velocities and the number of systems is too low for any statistically signihcant 
conclusions to be drawn. Although we should bear in mind that many uncertainties are 
involved and the current proper motion data do not allow any reliable conclusions to be 
drawn, we note that our results are consistent with the trend of close systems having higher 
center-of-mass velocities. 

We are grateful to Dimitrios Psaltis for invaluable discussions and comments on the 
manuscript. We thank Ron Webbink for a critical review of the manuscript, and also Fred 
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Fig. 1.— Post-supernova orbital separation, A, and eccentricity, e, as a function of the 
time, Tsn, elapsed since the supernova explosion, for PSR 1913-1-16 and PSR 1534-1-12. The 
vertical dotted lines he at Tsn = Tc, where Tc is the pulsar characteristic age (Table 1). 

Fig. 2.— Pre-supernova orbital separation, A^, mass of the exploding (helium) star, M^, 
and ratio, Ao/RHe,max, of the pre-SN orbital separation to the maximum radius reached by a 
helium star prior to its explosion, as a function of the time, Tsn, elapsed since the explosion. 
Quantities are plotted for PSR 1913-1-16 and PSR 1534-1-12 under the assumptions of circular 
pre-supernova orbits and symmetric supernovae. The vertical dotted lines lie at Tsn = Tc, 
where Tc is the pulsar characteristic age (Table 1). 

Fig. 3.— Ratio, Ao(l — eo)/RHe,max, of the pre-SN periastron distance to the maximum 
radius reached by a helium star prior to its explosion, as a function of the mass. Mo, of 
the helium star, for a complete range of pre-supernova eccentricities. Plots are shown for 
PSR 1913-1-16 and PSR 1534-1-12, under the assumption of symmetric supernovae explosions 
and Tsn = Tc. A minimum helium star mass for neutron star formation of 2 M 0 has been 
assumed. 

Fig. 4.— Limits on the pre-supernova orbital separation, Ao, and mass of the helium star. 
Mo of (a) PSR 1913-1-16 and (b) PSR 1534-1-12, for three different magnitudes of the kick 
velocity, 14. The vertical thick dotted line corresponds to the maximum orbital separation 
set by the geometrical constraint; the thick solid line corresponds to the minimum helium- 
star mass that can be accommodated in the orbit; thin lines correspond to the maximum 
possible pre-supernova helium-star mass (cos^ 0 = 0) and their position depends on the kick 
magnitude. Limits are calculated for circular pre-supernova orbits and for Tsn = Tc. 

Fig. 5.— Limits on the pre-supernova orbital separation, Ao, and mass of the helium star. 
Mo of (a) PSR 1913-1-16 and (b) PSR 1534-1-12, for three different values of the time, Tsn, 



elapsed since the supernova in units of the pulsar characteristic age, Tc- The thick solid 
line corresponds to the minimum helium-star mass that can be accommodated in the orbit; 
thin lines correspond to the limiting cases of cos^ 0 = 1 (vertical lines) and cos^ 0 = 0, and 
depend on the value of Tsn- Limits are calculated for circular pre-supernova orbits and for 
kick velocity magnitudes (a) 14 = 350km/s and (b) 14 = 300km/s. 

Fig. 6.— Limits on the pre-supernova orbital separation, Ao, and mass of the helium star. Mo 
of (a) PSR 1913-1-16 and (b) PSR 1534-1-12, for eccentric pre-supernova orbits, for Tg^ = Tc, 
and for different values of the kick magnitude, 14- 



max 



0 .2 .4 .6 .8 1 0 .2 .4 .6 .8 1 



0 .2 .4 .6 .8 1 0 .2 .4 .6 .8 1 

Tsn (lo" yr) Tsn (lo" yr) 















max 



0 .2 .4 .6 .8 1 0 .2 .4 .6 .8 1 



0 .2 .4 .6 .8 1 0 .2 .4 .6 .8 1 

Tsn (lo" yr) Tsn (lo" yr) 
























G 

s 





2.5 


3 


4 


4.5 

















TABLE 1 

Orbital and Pulsar Characteristics 



Units 

PSR 1913-1-16'= 

PSR 1534±12'' 

PSR 2303P46® 

PSR 1518±4904‘' 

PSR 2127±llCs 

Mtot" 

Mq 

2.82827 ±0.0004 

Binary Orbit 

2.68 ±0.003 

2.53 ±0.08 

2.62 ±0.07 

2.706 ±0.011 

Mp 

Mq 

1.442 ±0.003 

1.32 ±0.03 

~ 1.265'' 

~ 1.3^ 

~ 1.35'' 

M, 

Mq 

1.386 ±0.003 

1.36 ±0.03 

~ 1.265 

~ 1.31 

~ 1.35 

^cur 


0.617 

0.274 

0.658 

0.249 

0.681 

p 

cur 

sec 

27906.98 

36351.70 

1066136.65 

745978.02 

28968.36 

^cur 

R© 

2.80 

3.28 

30.5 

24.4 

2.83 

P ^ 

ms 

59.030 

Radio Pulsar 

37.904 

1066.372 

40.935 

30.529 

Ps 

10-'''s s-i 

8.626 

2.43 

569.3 

< 0.04 

4.99 

Tc 

yr 

1.1 X lO*' 

2.5 X lO*' 

3.0 X 10^ 

> 1.6 X 10'° 

9.7 X 10'' 

B 

10'°G 

2 

1 

80 

0.1 

1 

d 

kpc 

~ 5 

^0.5 

~ 2.3 

~ 0.70 

~ 9.8 


^Mtot- total mass, Mpi pulsar mass, Me- companion mass, Ccur: eccentricity, Peur: orbital period, Acur^ orbital semi-major 
axis 

^Ps spin period, Pg spin period derivative, Tc(= Ps/2Ps)'. characteristic age, B: magnetic field strength, d: distance from 
the Sun 

‘^Taylor & Weisberg 1982 

‘*Wolszczan 1991 
®Lyne & Bailes 1990 
*^Nice, Sayer, & Taylor 1996 
sPrince et al. 1991 

''For these systems, only the total masses have been derived. We estimate the individual masses by equating the masses of 
the two components. 







TABLE 2 

Visual Magnitudes of Possible Pulsar Companions'^ 




PSR 1913+16 

PSR 1534+12 

PSR 2303+46 

PSR 1518+4904 

PSR 2127+llC 

Main Sequence star 

~ 20 

~ 12 

~ 16 

~ 13 

~ 19 

White Dwarf 






Temperature (K) 

Age (yr) 






60000 

10® 

~ 26 

~ 18 

~ 22 

~ 19 

~ 25 

30000 

10" 

~ 27 

~ 19 

~ 23 

~ 20 

~ 26 

17000 

10® 

~ 29 

~ 21 

~ 25 

~ 22 

~ 28 

9000 

10® 

~ 30 

~ 22 

~ 26 

~ 23 

~ 29 


''We estimate visual magnitudes using the information about the companion masses, distances and positions on the sky. 
For main-sequence stars we adopt an F5 spectral type {My = 3.5), which is appropriate for the measured companion masses 
(Allen 1973), and for white dwarfs of different temperatures, and hence ages, we use the models by Bergeron, Wesemael, & 
Beauchamp (1995). We also apply a reddening correction using the maps of Burstein & Heiles (1982) and the extinction 
relation: Ay = E{B — V)/0.3. 


TABLE 3 

Pre-SN Parameters eor Symmetric Explosions 



PSR 1913+16 

PSR 1534+12 

PSR 2303+46 

PSR 1518+4904 

Ao" 

1.068 - 1.197Rq 

2.373 - 3.628 R© 

10.4 R© 

18 R© 


3.13-3.82 Mq 

2.1 - 2.3 M© 

2.9 M© 

2.0 M© 

■^o/ -^He,max 

0.12-0.36 

0.0084- 0.011 

0.31 

0.024 

^o(l 6o)/-RHe,max 

0.0024- 0.29 

0.0046 - 0.43 

0.022 - 2.5 

0.034- 3.16 


"Parameters are calculated for circular pre-SN orbits, over a range of DNS ages from 0 to IOtc (= Ps/2Ps). PSR 
2303+46 and PSR 1518+4904 experience no significant orbital evolution over these time scales. 

^’Parameters are calculated for a DNS age equal to Tc, over the entire range of allowed pre-SN eccentricities and 
He-star masses. 







TABLE 4 

Neutron Stars in Helium Common Envelopes 



PSR 1913-1-16 

PSR 1534-1-12 

PSR 2303-1-46 

PSR 1518-1-4904 

He-star Mass 

T 

^ ev 

3.2 - 2.9 Mq 

10 - 1 X 10” yr 

Evolution 

2.0Mq - 2.2 Mq 

9 X 10” yr 

3.2Mq - 2.9Mq 

1 X 10” yr 

2.0 Mq 

9 X 10” yr 



Accretion 



He-star Mass 

M 

TcoW 

4.2 - 2.5 Mq 

370 - 1.5 X 10” Meyr-i 
2.7 X 10-” -6.5 X 10-4 yr 

No Angular Momentum'^ 

2.0 M© - 2.25 M© 

0.9 - 250MQyr-’ 

1.1 -4.0 X 10-”yr 

2.5 M© 

2.5 M©yr-’ 

1.4 yr 

2.0 Mq 

0.11 MQyr-’ 

8.5 yr 

He-star Mass 

M 

Tcoll 

4.2 - 2.5 Mq 

160 - 200MQyr-4 

6.3 X 10“” - 5.0 X 10“” yr 

With Angular Momentum’’ 

2.0 Mq - 2.25 Mq 

0.01 - 75MQyr-’ 

100-0.013 yr 

2.5 Mq 
0.01M©yr-4 

110 yr 

2.0 Mq 
10-®MQyr-’ 

1.0 X 10® yr 


'‘We assume that the infalling material contains no angular momentum or is able to transport its angular momentum outward 
instantaneously. 

’’We assume that there is no angular momentum transport. We are thus limited to accretion close to the poles where angular 
momentum is insufficient to halt the inflow of the material. 







